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Abstract 

We present the UDF05 project, a HST Large Program of deep ACS (F606W, 
F775W, F850LP) and NICMOS (F110W, F160W) imaging of three fields, two of 
which coincide with the NICP1-4 NICMOS parallel observations of the Hubble 
Ultra Deep Field (HUDF). In this first paper we use the ACS data for the NICP12 
field, as well as the original HUDF ACS data, to measure the UV Luminosity 
Function (LF) of £ ~ 5 Lyman Break Galaxies (LBGs) down to very faint levels. 
Specifically, based on a V—i,i—z selection criterion, we identify a sample of 101 
and 133 candidate z ~ 5 galaxies down to ^50 = 28.5 and 29.25 magnitudes in the 
NICP12 field and in the HUDF, respectively. Using an extensive set of Monte 
Carlo simulations we derive corrections for observational biases and selection 
effects, and construct the rest-frame 1400 A LBG LF over the range Mmoo = 
[—22.2,-17.1], i.e. down to ~ 0.04 L* at £ = 5. We show that: (i) Different 
assumptions for the SED distribution of the LBG population, dust properties and 
intergalactic absorption result in a 25% variation in the number density of LBGs 
at ^ = 5: (ii) Under consistent assumptions for dust properties and intergalactic 
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absorption, the HUDF is about 30% under-dense in z = 5 LBGs relative to the 
NICP12 field, a variation which is well explained bjr cosmic variance; (iii) The 
faint-end slope of the LF is independent of the specific assumptions for the input 
physical parameters, and has a value of a ~ —1.6, similar to the faint-end slope of 
the LF that has been measured for LBGs at z = 3 and z = 6. Our study therefore 
supports no variation in the faint-end of the LBG LF over the whole redshift 
range z = 3 to z — 6. The comparison with theoretical predictions suggests that 
(a) the majority of the stars in the z — 5 LBG population are produced with 
a Top-Heavy IMF in merger-driven starbursts, and that (b) possibly, either the 
fraction of stellar mass produced in starburst, or the fraction of high mass stars 
in the bursts is increased towards the bright end of the LF. 

Subject headings: High Redshift Galaxies - Lyman-Break Galaxies - Reionization 
- Galaxy Formation - Galaxy Evolution - Cold Dark Matter 


1. Introduction 

Deep imaging with the Advanced Camera for Surveys (ACS) and the NICMOS camera 
onboard the Hubble Space Telescope (HST) has been instrumental in pushing the study of 
galaxy populations out to the reionization frontier of z ~ 6 and beyond. The development 
of the Lyman Break technique has been a milestone for the study of high redshift galaxies, 
as it allows to identify, from broad-band photometry alone, large numbers of star-forming 
but otherwise normal galaxies at early epochs: star forming galaxies at a given high redshift 
are identified from the drop in flux of, typically, 1-2 magnitudes, blueward of the detection 
passband, caused by the Lyman continuum break and the Lyman series blanketing induced 
by intergalactic neutral hydrogen clouds. While the conceptual idea was essentially worked 
out already 30 years ago (Meier 1976), it is only in the 1990’s that progress in instrumentation 
has made it possible to identify high-z galaxies using this technique (e.g. Guhathakurta et al. 
1990; Steidel & Hamilton 1992). Several spectroscopic follow-on studies of LBG candidates 
have proven the high efficiemry and reliability of the Lyman Break technique in identifying 
galaxies at early epochs (e.g. Steidel et al. 1999; Vanzella et al. 2006; Malhotra et al. 2005). 

Despite the availability of large samples of high redshift galaxies with low contamination 
from interlopers, the derivation of a Luminosity Function (LF), and particularly of its faint- 
end slope, remains difficult on account of the small volume probed by currently available 
ultra deep surveys such as the Hubble Ultra Deep Field (HUDF; Beckwith et al. 2006), and 
because of the many corrections which are often not well constrained. This is illustrated by 
the large scatter in the published LF parameters that have been derived from the same data 
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sets (e.g. Bunker et al. 2004; Yan & Windhorst 2004; Bouwens et al. 2006; Beckwith et al. 
2006). Mapping the evolution of the faint-end slope of the LF at high redshifts remains a 
major goal of astrophysical cosmology, as the faint-end slope is expected to be dramatically 
affected by reionization (e.g. Wvithe & Loeb 2006), and a solid detection of a change in slope 
would be a smoking gun for the transition between a partially and a fully reionized universe. 
As highlighted also by e.g., Beckwith et al. (2006), one important step in order, to improve 
the sensitivity to changes in the faint end slope is to derive LFs at different redshift in the 
most homogeneous way, applying similar techniques to similar data sets. 

This paper presents the UDF05 project, a 204-orbit HST Large Program of ultradeep 
ACS and NICMOS observations of multiple fields to study the evolution of the faint-end of 
the LF throughout the z ~ 4 — 8 redshift regime. The UDF05 was originally constructed so 
as to observe with ACS/WFC the two NICMOS parallel fields that were acquired while the 
HUDF was imaged with the ACS. These HUDF-NICMOS parallel fields (hereafter NICP12 
and NICP34, see Fig. 1), have a total exposure of ~ 2 x 10 5 s both in the F110W (broad 
.J) and F160W (H) filters,, and provide the deepest near-IR data available. Only relatively 
shallow ACS data were available at their location, and our additional UDF05 data for these 
fields was aimed at reaching a depth in the visible comparable to that of the NICMOS 
observations (see Table 1). Furthermore, the planned UDF05 pointings should have had 
an orientation such that, while acquiring the ACS data for the NICP12 and NICP34 fields, 
parallel NICMOS observations should have been taken for the original HUDF, so as to 
obtain a total of three fields with ultradeep imaging both in the visible and in the near-IR. 
Unfortunately, the transition of HST to the two-gyro mode and the first failure of ACS in 
June 2006 severely limited the number of available orbits at the required orientation, so that 
the NICMOS observations of the original ACS HUDF field ended up with only one third of 
the exposure time origina% planned. 

In our first utilization of the UDF05 dataset we study the faint-end of the LF of z = 5 
Veoe-dropout galaxies; specifically, we compare the z — 5 LF obtained from the ACS NICP12 
field with that derived from the original ACS HUDF, in order to assess the effects of cosmic 
variance on the z — 5 LF parameters. In Oesch et al. (2007, in preparation) and Stiavelli 
et al. (2007, in preparation) we will extend our exploration of the realistic uncertainties in 
the faint-end slope of the LF of Lyman Break Galaxies to the z = 6 and z W 6 regimes, 
respectively, probing well into the expected reionization epoch and thus constraining its 
impact on galaxy formation. 

This paper is structured as follows. After the description of the data, corrections for 
noise correlation and source detection (§2), we describe our selection criteria (§3), and the 
simulations that were performed to estimate the completeness of our source catalogs (§4). 
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The observed surface densities are derived by including the effects of photometric errors (§5) 
and the LF is computed under different assumptions about the underlying SED distribution 
of the LBG population as well as different prescriptions for intergalactic hydrogen absorption 
(§6). In section §7 we present a theoretical estimate for the effects of cosmic variance on 
our results. In section §8 we discuss the evolution of the faint-end of the LF of star forming 
galaxies in the 2 = 3 — 6 redshift regime, and compare our observed z — 5 LF with theoret- 
ical predictions based on a standard 1 A Cold Dark Matter (A-CDM) universe. Section §9 
summarizes our main results and conclusions. A comparison of our LF measurements with 
published observational estimates is shown in Appendix A. 


2. Data 

The ACS WFC data on the HUDF are the same as described by Beckwith et al. (2006) 
and include exposures in the F435W (hereafter £? 435 ), F606W (hereafter Vfog), F775W 
(hereafter *775)5 and F850LP (hereafter z 850 ) filters for a total of 400 orbits. We used the 
publicly released images without further processing. 

The UDF05 ACS WFC data for the NICP12 field in the F606W, F775W, and F850LP 
filters amounted to a total of 102 orbits; Table 1 lists the specifics of these observations. The 
center of the NICP12 field is located at RA(2000)=03:33:03.60, Dec(2000)=-27:41:01.80; 
Figure 1 show's all UDF05 pointings relative to the HUDF ACS and NICMOS-parallel data. 
The UDF05 ACS observations of the NICP12 field employ a larger dithering step than those 
adopted for the original HUDF ; this choice was made in order to obtain an improved sub- 
pixel sampling for the NICMOS parallel data that were acquired simultaneously, relative to 
the original HUDF NICMOS-parallel observations. 

The exposure times for the UDF05 observations w-ere designed so as to achieve a constant 
AB-magnitude depth in all filters for all fields. The individual 1200s exposures were combined 
and re-pixeled to a scale of 30 mas using the task ’’multidrizzle”; the task ’’tweakshifts” 
was used to optimize the shifts between different images. Superdarks and superbiases were 
obtained combining several hundred single frames to increase the S/N. Finally, a herringbone 
artifact was removed from the data by processing the images with a Fourier filter, and the 
electronic ghost present in the ACS images was eliminated by means of an empirical model. 
The model is built upon the fact that the electronic ghost of a source is located in quadrants 
different from the one containing the source and on the fact that ACS/WFC is read out from 


1 Throughout this paper we adopt the concordance cosmology: Qj{ = 0.3, Qa — 0.7, Hq = 

kms -1 Mpc -1 , i.e. h = 0.7. Magnitudes are given in the AB system (Oke & Gunn 1983). 


70 
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each corner so that if one dithers a source closer to its readout corner the ghost will move 
in the opposite direction closer to its own corner. Thus, when we drizzle together a dithered 
pair of images we align the sources and smooth out the ghosts. Our algorithm (conceived by 
L.E. Bergeron) is based on the idea of flipping each quadrant before drizzling. This simple 
step coadds the ghosts and smooths out the sources and when applied iteratively with proper 
masking can be used to obtain clean images of the ghosts which are then subtracted from the 
images. These last two steps were developed after the first public release of the UDF05 data 
(version 1 , ”vl”). The newly processed images of the UDF05 that we use in our analysis 
represent an improved ’’version 2 ” (v 2 ) of the final products. 

The PSF FWHM of the UDF05 final images are about O'.' 1 in all passbands. In the 
UDF05 NICP12 images, the 10 a limiting magnitudes in apertures of 0/15 radius are 28.49, 
28.44, and 28.47 respectively in l/ 60 6 , 2775 , and z$ so (see Table 1), i.e. the planned uniform 
sensitivity in all passbands is achieved. For the HUDF the B. 435 , 1-606 ^ *775 5 and zgo 0 10 a 
limiting magnitudes in apertures of O'/I 5 radius are 29.38, 29.81, 29.43, 28.73 respectively. 


2 . 1 . Noise Correlation 

A point kernel was used for all passbands when drizzling the data, with the exception of 
the NICP12 Vgoe image, for which a non-zero pixfrac had to be used, as the dithering pattern 
was inadequate to properly sample the pixels, i.e. for this band, the pixels were not shrunken 
to a point before drizzling to the output frame and, as a result, this image is expected to 
show a significant noise correlation. However, the noise level measured by computing the 
rms of cleaned, block averaged areas of the images, or bj^ computing the total flux in 
the auto-correlation peak shows a higher degree of correlation than the one due to dithering 
which is well characterized by the model published by Casertano et al. (2000). Moreover, 
some noise correlation is also measured in the i - 75 and 2:350 images which were not expected 
to have any. The unexpected correlated noise seen in these images is of the same amplitude 
as the additional component seen in the f / 60 6 image and we believe that this residual noise 
correlation is introduced by the reference files (flat, fields, biases, darks) and, possibly, by 
unresolved background sources. Any feature present in a given pixel in the reference files 
is propagated to a larger area on the final science images by the adopted dithering pattern. 
Neglecting this noise component would lead to overestimating the depth achieved by the 
observations; for example, in the NICP12 field, the noise would be underestimated by 30% 
for galaxies with an area of 0.3 x 0.3 arcsec 2 . 

Relative to our final UDF05 images, the released ACS HUDF images that we use in this 
paper are more severely affected by correlated noise, due to the presence of the uncorrected 
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herringbone noise and electronic ghost in ACS (as these effects had not been characterized 
at the time). This explains why, despite the factor of about two in exposure time, the HUDF 
zgso image is only 0.26 mag deeper than our NICP12 z 850 image. 

The rms maps produced by multidrizzle were rescaled in order to match the expected 
flux errors of galaxies with a measured area of 80 pixels, which corresponds to the median 
size of our dropout galaxies. No simple renormalization is optimal for all source sizes. The 
adopted choice underestimates the S/N for point sources and overestimates it for sources 
that are more extended than the considered area; it provides however on average an optimal 
renormalization for the tj’-pical galaxies that we study in this paper. Simulations are used to 
assess proper errors as a function of the sources’ physical parameters. 


2.2. Source Detection 

Source detection and photometry measurements were done with the software tool SEx- 
tractor (Bertin & Arnouts 1996), which was run in double image mode with 2 r 85 o as the 
” detection band”. Despite the fact that the i r75 images for both the HUDF and the NICP12 
data reached a similar if not deeper magnitude limit than the z 850 band, the la.tt.er was 
preferred for the source detection because of the larger uncertainties in an ? T75 -selected LF, 
which are driven by the redshifting of the Lycc forest into the i 7T5 passband at z > 4.8, 
significantly decreasing the observed fluxes. 

The SExtractor parameters were optimized so as to maximize the number of detected 
galaxies while minimizing the number of spurious sources. For NICP12, the detection thresh- 
old was set to be a minimum of 10 connected pixels 0.55 a above the background. This led 
to 5,734 objects within the central, highest S/N 11.9 arcmin 2 of the z 850 NICP12 image. 
Similar SExtractor configurations returned 6,336 objects in the central 11.2 arcmin 2 of the 
HUDF z% so image. 

In order to estimate the contamination of spurious detections in the source catalogs, we 
cleaned the images from the detected sources, and reran SExtractor using the same configu- 
ration parameters on the ’’inverted” (i.e. , multiplied by -1) images. This contamination was 
found to be < 0.5%. A final visual inspection was performed to clean the catalogs of residual 
spurious sources and stellar diffraction spikes. The source magnitudes were corrected for the 
small amount of Galactic extinction expected in the field, E(B — V) = 0.007 and 0.008 in 
NICP12 and HUDF, respectively (Schlegel et al. 1998). 

Colour measurements were derived from the galaxy ISO magnitudes (apertures matched 
to the detection isophotes), since these are closely matched to the higher S/N parts of the 
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objects and, as shown by our simulations (Section 4), include much less noise than the 
SExtractor AUTO magnitudes which we use for the total magnitude of a source in the 
detection band. This is especially important for the very compact and faint, sources that we 
study in this paper. 


3. The Selection of z ~ 5 Lyman-Break Galaxies 
3.1. The Color-Color Selection Criterion 

The redshift evolution of synthetic galaxy models, properly attenuated by intergalactic 
hydrogen absorption, can be used to identify regions in a color-color diagram which efficiently 
disentangle galaxies in a specific high redshift window from lower-redshift interlopers. 

For the identification of a V’eos-dropout sample (see Figure 2), this selection window has 
been optimized to be (see Giavalisco et al. 2004; Beckwith et al. 2006): 


Veoe - ir~5 > min [ 2 , 1.5 + 0.9 • (i r7 5 - z 850 )) g§) 

he06 — *775 >1-2 (2) 

*775 ~ '*850 < 1-3. §§J 

Galaxies in the range z ~ 4.5 — 5.7 are well identified by this color-color criterion; the 
reliability of the method has been shown to be very efficient (^ 90%) by spectroscopic 
follow-up surveys (Vanzella et al. 2006). 

The elimination of low redshift interlopers is facilitated by the availability of the B 4:i 5 
passband for the HUDF sample. This is however not available for the NICP12 data. We 
therefore used the HUDF data to check the impact of including or excluding the additional 
S /N ( B 435 ) < 3 criterion that was adopted by Giavalisco et al. (2004); Beckwith et al. (2006) 
in order to minimize the contamination from low- z interlopers. Only one object was rejected 
from our HUDF dropouts catalog on the basis of this 1 % 135 -band constraint. Therefore, for 
uniformity with the analysis of the NICP12 field, for which the B 435 is not available, we did 
not include any B 4 35 -band constraint in our selection of the HUDF Veoe -dropouts sample. 

The stellar library of Pickles (1998) shows that stars are expected to lie within a very 
well defined region in the Veoe — *775, *775 — ~sso diagram. Indeed, the few obvious stars in 
our field lie within this area. To further optimize the removal of stars from our catalogs, 
the SExtractor stellarity parameter was also considered. This parameter is determined by 
application of a neural network approach to estimate the likelihood that an object is a 



point source. However, this method is not very reliable at faint magnitudes, and therefore 
cannot be used to identify stars close to the magnitude limits. Inspection of the Veoe — £775 
vs. ?~75 — zg 5 o diagram and of the half-light radius (rg/ 2 ) vs. z $ 50 diagram showed that 
sources with stellarity > 0.85, z 850 < 27 and 7 * 1/2 < O'/I are most likely stars in the NICP12; 
lienee all such objects were excluded from our analysis. A similar criterion was applied to 
the HUDF for z 8 5 o < 28. Only two and one of such sources were within our I/fog-dropout 
selection window for the NICP12 and HUDF samples 2 , respectively; these objects however 
had already been visually identified as stars. Although unlikely, stars fainter than z 8 50 ~ 27 
and 28 may remain a (small) source of contamination in our NICP12 and HUDF catalogs, 
respectively. 

In Fig. 3 the evolutionary tracks of different galaxy types are shown in the V m — i 750 
versus f-75 — z 850 color-color diagram. The tracks are built using the population synthesis 
models of Bruzual & Chariot (2003). The attenuation by intergalactic hydrogen causes es- 
sentially all the galaxy types at redshifts above WL5 to enter the selection window. Elliptical 
galaxies at z ~ 1.5 are a source of contamination in the selection criterion due to their red 
colors; this contamination is however expected to be small in our V' 60 6-dropout samples, given 
the expected low surface densities of these galaxies. In fact, none of the seven spectroscopi- 
cally identified passively evolving galaxies at a > 1.4 in the HUDF from Daddi et al. (2006) 
were selected as Vjoerdropouts. 


3.2. Contamination from Low-Redshift Interlopers 

We have estimated the possible fraction of interlopers by applying our selection criterion 
of Eqs. (l)-(3f to a library of ~ 3000 synthetic SEDs built upon Bruzual-Charlot (2003) 
models, adopting the LF derived by Steidel et al. (1999) at z = 3 and no evolution. The 
models include the effects of intergalactic absorption (Madau 1995), and span a wide range 
of metallicities (0.04 — 2.5Z©), dust reddening, emission lines, and different star formation 
histories. These included simple stellar populations, models with continuous star formation 
at constant metallicity, models with self-consistent enrichment during the continuous star 
formation (estimated with a closed box model and with a model with infall of zero- metallicity 
gas), and two-burst models combining with five different mass ratios an old (> 50 Myr) 
stellar population with a young (1 — 10 Myr) star forming population that includes self- 
consistent nebular line and continuum emission. For each model, Aversions with and without 


2 All other sources in the HUDF stellar catalog of Pirzkal et al. (2005) that we do not identify as stars lie 
outside and very far from the border of our color-color selection window' . 
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dust reddening were computed; eight values of reddening were implemented assuming a 
galactic dust screen with Ay logarithmically spaced between 0.05 and 6.4. Emission lines 
and nebular continuum were computed in the approximation of case-B recombination for 
Hydrogen and Helium, based on the ionizing flux from the stellar SED; Oxygen and Carbon 
abundances were derived from an analytical interpolation of a grid of calculations performed 
with the code CLOUDY (Ferland et al. 1998). Gas temperatures were derived as a function 
of metallicity. 

The resulting redshift distribution is shown in Fig. 4. The non-uniform distribution 
of the histogram is due to SEDs that satisfy the selection criterion onty for portions of the 
redshift interval. For instance a 1 Gyr old solar metallicity model with modest extinction 
satisfies the criterion only for a > 5.5. We identify three possible classes of low redshift 
interlopers at z ~ 0.6, ~ 1.6, and ~ 3.2. They are generally very reddened SEDs with either 
emission lines or an underlying old population. The total interloper fraction from these 
objects is found to be 23% (17% if one considered only z < 4 as interlopers). This is only a 
rough estimate, as it does not account for the different galaxy populations having different 
and possibly evolving LFs. However, previous works show that the utilization of this library 
overpredicts the number of interlopers (e.g. with respect to Malhotra et al. 2005); therefore, 
the real fraction of interlopers is possibly smaller than our conservative estimate. 


4. Quantification of the Catalogs’ Completeness: Simulations 

The catalogs produced by SExtractor are subject to certain biases, which depend 
strongly on the chosen parameters. In addition to the detection settings, i.e. , number of 
connected pixels and sigma-threshold above the background, the parameters that control 
the deblending and cleaning of sources are of particular importance. The correct interpreta- 
tion of the LF that is obtained from the S Extractor-based catalogs requires a quantitative 
assessment of these biases. To this end, we performed two independent sets of simulations, 
the first based on idealized galaxy profiles and the second on the observed galaxies them- 
selves. The simulations are described in detail in Appendix B. Briefly: 

1. In the first set of simulations (which we refer to as the ” idealized-galaxies” simulations), 
we constructed an extensive set of exponential- and deVaucouleur-profile galaxies, that 
were appropriately PSF-convolved and inserted at random positions into the NICPI2 
and HUDF images. SExtractor was then re-run with identical parameters as for the 
compilation of the original data, catalogs, and the resulting object catalogs were used 
to identify detected, undetected or blended sources. 
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2. In the second set of simulations (which we refer to as the ’’dimmed-galaxies” simula- 
tions), we generated a set of test images by co-adding to each original image a copy of 
itself, shifted by a few pixels and dimmed by a given amount. The composite images 
were rescaled to the S/N of the original image, and SExtractor was run on these final 
test images using the same parameters as for the extraction of the original data cata- 
logs. This second set of simulations provides a complementary analysis to the above, 
as it offers the advantage of working with real rather than idealized galaxies, and does 
not depend on uncertainties in the adopted PSF. 


4.1. Input versus Output Properties 

A problem in establishing the completeness of the source catalogs is that the true mag- 
nitudes of the sources are expected to be typically brighter than the measured ones, due to 
losses of flux from the wings of the galaxy light distribution. Sirianni et al. (2005) measured 
the flux losses for point-sources in the ACS camera filterbands and found that they strongly . 
depend on the SED of the source, due to a stronger scattering halo of the PSF at longer 
wavelengths. Those authors find offsets of ~ 0.1 mag for and i 77 5 , respectively, and 
~ 0.15 mag for z 85 o in apertures of radius 0.4". To take into account the extended nature of 
our catalog objects, we used our two sets of simulations to characterize the relation between 
input (true) parameters and output (measured) quantities. 

Idealized deVaucouleur- and exponential-profile galaxies were generated with random 
,2850 magnitudes between 24 to 31. Two separate sets of tests were performed, the first 
using sizes uniformly distributed in the range 0.05 to 0.5 arcseconds, and the second using 
sizes distributed according to a. lognormal distribution centered at 0.25 arcseconds with 
a = 0.3. The results were essentially independent of the adopted input size distribution 
and in the follouung we only use the simulations with uniformly distributed input sizes. 
Measuring the source magnitudes within an elliptical aperture of 2.5 Kron radii (SExtractor 
AUTO mag), compact sources were detected with 50% probability down to 28.75 and 29 
magnitudes in the NICP12 field and in the HUDF, respectively. As expected, the detected 
fluxes were typically smaller than the input ones. The theoretically expected flux losses for 
the deVaucouleur and exponential profiles are 9.6% and 4.0%, respectively, corresponding to 
0.11 and 0.04 magnitude offsets. Since, however, the Kron radii measured by SExtractor are 
also underestimated, the offsets are expected to be even larger than those above. Indeed, in 
our simulations we found a magnitude-dependence for the amount of dimming, from ~ 0.3 
mag and ~ 0.07 mag at % 50 = 24.25 to ~ 0.6 and ~ 0.2 mag at z 850 = 28.25 for the two 
kinds of profiles, respectively. 
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We also used the dimmed-galaxies simulations to obtain an independent estimate. In 
this case, the differences between dimmed and original magnitude are expected to be smaller 
than in the case of idealized galaxies. The offset was indeed marginally smaller than for the 
exponential-profile idealized galaxies, i.e., in the range Amag ~ 0.04 — 0.12. 

For very compact galaxies (r i / 2 < 0715), the measured sizes were generally overpredicted 
due to PSF blurring, while the opposite was true for larger objects. These biases showed a 
strong dependence on magnitude. The faintest idealized deVaucouleur galaxies were typically 
measured to have ri /2 ~ 071 (i.e. to be unresolved at the resolution of the HST), independent 
of their original size. This offers a challenge for studying the size evolution of LBGs with 
redshift. 

Furthermore, since we measure the colors bjr relying on the SExtractor dual-image mode, 
i.e. using the z 850 detection apertures, it is not adequate to apply a constant offset to the 
y 606 and £775 magnitudes. The color measurements depend in fact on the z&o flux as well as 
on the intrinsic color. Both the £775 — z s 50 and the Veo 6 — *775 colors are affected in such a 
way that red sources have even redder, and blue sources even bluer, measured colors. This 
is not including the dependence of the PSF halo on the color. 

In the light of these complex and uncertain corrections for dimming effects, we chose 
not to implement any correction to the measured catalog magnitudes. Thus, we note that 
our magnitude measurements might be underestimated by an amount, as estimated from 
our tests, of up to about 0.2 magnitudes. 


4.2. Detection Completeness 

There are two possible reasons why galaxies could be missing from our catalogs. First, 
a galaxy might not be detected due to a too low surface brightness; second, it could get 
blended together with another source and produce a false catalog entry. These two effects 
are strongly dependent on the SExtractor parameters. 

We used both the idealized-galaxies and dimmed-galaxies simulations to estimate the 
strength of these effects in our catalogs and final samples of y 6 06 -dropout galaxies. Specifi- 
cally: (a) We associated idealized-galaxies with similar output (SExtractor) magnitudes and 
sizes to real galaxies randomly selected from our original catalog. We then estimated the 
completeness corrections for each magnitude bin by comparing the number of sources in each 
bin of SExtractor magnitude with the total number of idealized-galaxies of different input 
(i.e., theoretical) magnitudes and sizes that contributed to that specific bin of output mag- 
nitude. (b) Using the dimmed-galaxies simulations, we computed the magnitude-dependent 
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completeness corrections by re-scaling the number of galaxies in each bin of output mag- 
nitude (i.e., recovered for the dimmed galaxies by SExtractor) to the number of sources of 
different input magnitudes (i.e., the theoretical values of the dimmed magnitudes) that con- 
tributed to that specific bin of output magnitude. Figure 5 shows the results of our tests. For 
both fields, the two independent estimates for the completeness in our samples, derived from 
the idealized- and dimmed-galaxies simulations, are in good agreement. At the magnitude 
limits of our selections, the NICP12 and HUDF catalogs are expected to be about 65% and 
50% complete, respectively. 

The probability that a given source is blended by SExtractor with another object de- 
pends on the magnitude of the source. In our simulation tests, when the SExtractor de- 
tection blended a simulated source with a real galaxy, we considered the simulated object 
as a detection when this was brighter than the original galaxy. Figure 5 also includes the 
incompleteness due to source blending; this is the reason why even our brightest sources are 
only complete at the 95% level. 


5. Number of Vgo6-dropout Sources 

The application of the selection color window discussed above to our NICP12 and HUDF 
catalogs down to 28.5 and 29.25 magnitudes (with a S/N > 6 cut) produces a total of 101 and 
133 candidate z ~ 5 LBGs, respectively; Note that for sources whose V m ISO fluxes were 
fainter than their 2a-errors, we replaced the measured Vsoe fluxes with the corresponding 20- 
values, in order to get a reliable lower limit on their V WQ — £775 color. 

A bootstrapping technique was applied to estimate the impact of photometric scatter 
and the effect of the binning in 0.5 mag bins on the selected number of Vgog-dropout sources; 
objects in our catalog were assigned photometric errors derived from the simulations accord- 
ing to their total magnitudes and ISO colors, before the selection of Vgog-dropout sources 
was repeated. 

The idealized- galaxies simulations were used to estimate the fraction of interlopers f int} 
i.e. the fraction of sources that were included in z = 5 LBG sample with an expected 
intrinsic color that does not match the selection; the f int value was used to correct the 
detected number of LBGs so as to take into account this contamination by interlopers. The 
detected and corrected numbers of LBGs in the different bins of magnitudes for both the 
NICP12 and the HUDF are listed in Table 2, together with the applied corrections. 

Tests were performed to ensure that the final number of sources in the LBG catalogs were 
independent of the adoption of the 2a U 306 flux limits for faint Vgoe sources. In particular, the 
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calculations were repeated adopting la and 3 a limits for the flux substitutions; no significant 
differences were found in the final LFs, as the probability of selecting a galaxy at any given 
redshift changed accordingly (see next section). 


6. The Luminosity Function of z — 5 Lyman-Break Galaxies 

The UV-continuum LF was derived from a maximum likelihood fit of a Schechter (1976) 
function to the observed numbers of F 606 -dropout galaxies in bins of apparent magnitudes, 
JVj, with: 


rm.-xt 

Ni~ / dm N(m;<l>*,M*,(x) = (4) 

J mi 

f°° dV f mu 

/ dz—— / dmp(m,z)<i){M{'m,z)\<j)*.M*,Ci) 

Jo dz J mi 

Here nq and m u are the lower and upper bounds of the i - th magnitude bin, centred 
around m*. The factor p(m, z) is the probability that a galaxy with magnitude m at redshift 
z is detected and selected to be az = 5 LBG. 

An ’’effective volume” (14ff) technique was also used to compute the stepwise LF in bins 
of absolute magnitude under the assumption that M (m, z) is slowly-varying with redshift, 
i.e.: 

f°° dV 

Ni ~ z)) / dz—p{m h z ) (5) 

= <j>(M(m. h z)) V; ff (m,;) (6) 

Over the redshift range of interest, z = 4.5 — 5.7, the absolute magnitude at any given 
apparent brightness changes by about 0.5 magnitudes, i.e., an amount comparable with the 
size of our magnitude bins; the resulting LF is thus generally less robust, 3 than the one 
computed using Eq. 4. 

In order to properly constrain the value of M» in our best fit Schechter function, we adopt 
this version of our LF for matching it to the Yoshida et, al. (2006) data which constrained 
the z = 5 LF at bright magnitudes in the Subaru Deep Field (SDF) (see Section 6.2). 


3 Not,e that the situation would be exacerbated for an «77s-selected (rather than zgso -selected) Veoc;- 
dropout sample, since the differences in absolute magnitude over the same redshift range could be up to ~2 
magnitudes, i.e., much larger than the size of the bins, due to intergalactic Lya- absorption. 
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The K-corrections required to transform the observed magnitudes into rest- frame 1400 
A absolute magnitudes, were computed assuming a 200 Myr old, continuously 

starforming galaxy template, reddened by E(B — V ) = 0.15 with a Calzetti et al. (2000) dust 
extinction relation (hereafter referred to as ’’starburst” dust extinction). The differences in 
the K-corrections for alternative reasonable assumptions on the stellar population properties, 
including significantly younger stellar populations, -were checked and found to be negligible 
(< 0.05 mag). 

The probability that a galaxy with magnitude m at redshift z is detected and selected 
to be a z = 5 LEG can be decomposed into the product of two components, i.e. , the 
completeness-probability C(m) and the selection-probability S(m,z ), i.e. p(m,z ) = C(m ) • 
S(m, z). We discuss this in detail below. 


6.1. The Selection Probability S(m,z ) 

The selection-probability as a function of magnitude and redshift was computed includ- 
ing the scatter in the photometry, which was derived from our simulations, and replacing with 
the 2 — a flux value the V 60(i magnitudes fainter than this limit. The input assumptions for 
the distribution of SEDs in the Vg 06 -dropout population were based on the findings of Steidel 
et al. (1999). These authors inferred a distribution of spectral slopes from their sample of 
spectroscopically-confirmed z ~ 3 LBGs by assuming a 1 Gyr old, continuously starforming 
galaxy SED with varying dust content according to the starburst dust extinction relation 
(see also Adelberger & Steidel 2000). The E(B — V) distribution -was found to be roughly- 
approximated by a Gaussian with a median of 0.13 and a a — 0.12. Note that this results 
in a small fraction of the galaxy population which is fitted by non-physical, negative values 
of E(B — V); this might indicate the presence of strong Lycc emission for some galaxies. 
Nevertheless, this ’’effective” distribution of E(B — V) is in good agreement with the typical 
dust attenuation values found by Papovich et al. (2001), who used synthetic SED fitting to 
z ~ 3 LBGs to infer the properties of these galaxies. Since an age of 1 Gyr might be too old 
for z ~ 5 galaxies, we used similar models with a 200 Myr old, constantly starforming SED 
(Bruzual & Chariot 2003). The selection probability was found to be virtually independent 
of the assumed stellar age; this is understandable in the light of the dominant contribution 
from the youngest stellar populations to the UV emission. 

To explore the effects of different assumptions for the intergalactic hydrogen absorption, 
the SEDs were attenuated according to the calculations of both Madau (1995) and Meiksin 
(2006). In the latter case, the selection probabilities are slightly reduced, as a consequence 
of the smaller attenuation and thus smaller reddening of the SEDs. 
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We estimated the impact on the selection probability of assuming that there are no 
Lyman-Limit systems along the line of sight attenuating the observed Vqq 6 flux, a situation 
which is expected to happen with a probability of only 6% up to z = 5. This effect was 
found to be negligible, i.e., the selection probability was reduced by less than 5% around 
z ~ 5.2. In the remainder of our analysis we therefore assume the estimated mean values for 
the intergalactic absorption. 

We also tested alternative distributions for the UV slopes. Verma et al. (2007) used 
stellar population synthesis models to derive ages, masses, and dust attenuation properties of 
z ~ 5 LBGs. Using the Small Magellanic Cloud (SMC) dust extinction law, both Papovich et 
al. (2001) and Verma et al. (2007) find that the typical UV attenuation is smaller than in the 
case of the starburst extinction law. The exact values somewhat depend on the metallicities 
and Initial Mass Functions (IMFs). To estimate the selection probability, we adopted a 100 
Mvr old, continuously starforming template reddened by a Gaussian distribution of E(B—V ) 
with mean= 0.05 and a = 0.04, which is approximately in agreement with the distributions 
found by Papovich et al. and Verma and collaborators. Also in this case we explored the 
effects of changing the prescriptions for the hydrogen intergalactic absorption. Differences 
of order 10% or less were found with respect to the previous estimates. 

The selection probabilities for the NICP12 field, as derived under the various different 
assumptions, are shown in Fig. 6 for a bright and a faint magnitude bin. In general, 
since the intergalactic absorption relation derived by Meiksin (2006) is smaller than the 
corresponding calculation of Madau (1995), galaxies are less reddened assuming the former 
prescription; consequently, at any redshift of interest, they lie closer to the boundaries of our 
color selection. This leads to a reduced probability that they are selected as LBG candidates. 
The effect is exacerbated when the SMC extinction law is assumed: a considerable drop in 
the redshift distribution is expected around z ~ 5.3, and the lower limit for the selection- 
redshift is slightly increased. Note that the adoption of the 2 — a lower limit for the Vgoe fluxes 
leads to a considerably different redshift window for the fainter magnitude bin: the reddest, 
highest-redshift objects are not detected. Furthermore, due to the larger photometric scatter 
at fainter magnitudes, the redshift selection smears out to substantially lower redshifts, down 
to z ~ 4. Note that we have not applied any correction for incompleteness in these plots; 
therefore, the drop in selection probability at fainter magnitudes that is observable in the 
Figure is exclusive^ due to the errors in the color measurements. Similar results hold for 
the HUDF, for which however the selection probabilities are somewhat larger, due to the 
smaller uncertainties in the color measurements (see bottom panel of Fig. 6). 
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6.2. The Impact of Different Physical Assumptions on the LF 

As we discuss below, different assumptions for physical input parameters such as the 
dust extinction law and the hydrogen intergalactic absorption relation have a non-negligible 
impact on key quantities such as the effective volume and thus the average luminosity den- 
sity. It is thus of paramount importance that LFs derived from different data sets are first 
rescaled to identical assumptions for SED distributions, dust properties and intergalactic 
extinction, before they can be compared to study any possible evolutionary effect in the 
galaxy populations. In our analysis we adopt the Madau (1995) and starburst relations as 
our ” fiducial” assumptions for the hydrogen absorption and dust extinction, respectively, as 
these are widely used in most previous works and allow for a direct comparison with the 
published results. We explore however below the effects on the LF parameters of varying 
this assumptions. 

Our pencil-beam observations do not properly constrain the value of M* in the LF: 
given the rarity of very bright objects, the small volumes probed by the NICP12 field and 
the HUDF are not sufficient to adequately sample the knee of the Schechter function. We 
therefore rely on the SDF z = 5 LF of Yoshida et al. (2006; see also Appendix A.l) to 
constrain the bright-end of our U e ff LF. These authors used similar fiducial assumptions for 
the intergalactic absorption and dust extinction relations, and similar techniques to correct 
for the observational biases, to those that we adopt in this paper. 

The SDF LF is based on the Ri'zf SUBARU filters, and is sensitive to a redshift window 
that is only slightly narrower than that of our data; it reaches down to z! = 26 mag, i.e., ~ 3 
magnitudes brighter than our measurements. In combining the SDF and NICP12/HUDF 
LFs, we neglected small K-corrections from 1500 A to 1400 A; the V e g SDF LF nicely matches 
our Left NICP12 LF in the region of overlap. The total Vgir LF, that merges the SDF estimate 
with our NICP12 and HUDF LFs, spans over the 5-magnitudes range Mu oo — [ — 22.2, —17.1] 
(Figure 7). 

Adopting the completeness corrections derived from the idealized-galaxies simulations, 
the lowest % 2 , best fit of a Schechter function to the total I4ff LF provides the following 
parameters: </>* = 0.9loJ • 10~ 3 Mpc -3 mag -1 , M* = — 20.78tg;J® mag, a = —1.54^0;];° (see 
Figure 7 and Table 4). As expected, the value of M, * is largely influenced by the SDF 
LF, while the slope a is constrained by our fainter data. The adoption of the completeness 
corrections derived from the dimmed-galaxies simulations, and of different SED distributions, 
dust extinction laws and intergalactic absorption prescriptions, has no significant impact on 
this result. 

In particular, all tests performed by varying these quantities, independent of whether 
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they include one or both of our two deep fields, return a value of M* that, within the errors, 
is consistent with M» = —20.7. We therefore fixed M* to this value in order to highlight the 
effects of varying the input physical parameters on the LF computed according to Eq. (4) 
for the HUDF and NICP12 only. 

The results obtained with fixing the value of M* = —20.7 and not using the SDF for 
the fits are listed in Table 3. Only marginal differences are generally observed: within 
the errors, the derived parameters are consistent with each other. Note however that the 
average luminosity density increases by about 20-30% when the Vgoe-dropout population is 
described by the SMC dust extinction law and the Meiksin intergalactic absorption relation, 
relative to the more widely used combination of prescriptions which we have adopted as our 
fiducial model. This is due to a change of a similar order of magnitude in the resulting 
effective volumes. We note that the uncertainties in the completeness, estimated from our 
two complementary sets of simulations, have a comparably small effect on the faint-end slope 
of the LF as the different physical assumptions. 


7. The Role of Cosmic Variance 
7.1. The Underdensity of the HUDF 

On average, independent of the assumed input parameters, the LF derived from the 
HUDF results in a lower average luminosity density: The HUDF appears to be "under- 
dense” with respect to the NICP12 field (see also Fig. 7). By integrating the LF from 
Table 3 over the magnitude range of overlap between the HUDF and the NICP12 field, 
i.e. , M 1400 = [—21, —18], the amount of this underdensity is about 30%, independent of the 
specific input physical assumptions. Although a larger number of low-redshift interlopers in 
the NICP12 sample cannot be ruled out, (a) the homogeneity in the procedures with which 
both data sets have been handled and analyzed, and (b) theoretical expectations concerning 
the amount of scatter amongst volumes of order of those sampled with the HUDF and 
NICP12 field, both support the interpretation that the under-densitv in z ~ 5 LBG of the 
HUDF has a physical origin. Note that, in their analysis of the f-dropouts population, 
Bouwens et al. (2006) find the HUDF under-dense at the bright end, also by about ~ 30%), 
relative to the GOODS field. Thus, it appears that the ~ 30% under-density of the HUDF 
extends at least over the volume encompassing the entire z = 4.5 — 6.5 epoch. 
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7.2. Quantification of Cosmic Variance Effects in Our Study 

In order to estimate the impact of cosmic variance on our measurements, we use a 
5 12 3 -p articles, dark-matter-only cosmological simulation in a box of 100 Mpc/h side, that 
results in a particle mass of 5 • 10 8 M e /h. Simple prescriptions to populate the dark matter 
halos with luminous galaxies, and a realistic beam-tracing algorithm were used to study the 
field-to-field variations in the number counts of star forming galaxies at the epoch of interest. 

The numerical simulation was carried out with the public PM- Tree code Gadget-2 
(Springel 2005) adopting a ug = 0.75 value. The initial conditions were generated using a 
code based on the Grafic algorithm (Bertsehinger 2001) with a, AC DM transfer function com- 
puted according to the fit of Eisenstein & Hu (1999), and assuming a spectral index n = 1. 
Dark matter halos were identified in the simulations snapshots using the HOP halo finder 
(Eisenstein & Hut 1998); the redshift interval 4.5 — 6 was sampled in steps of A z — 0.125. 
The simulation reproduces a volume about 73 times larger than the effective volume probed 
by each of our two fields for the V-dropouts population (« 5.7 x 5.7 x 420 (Mpc/h) 3 ). 

Source catalogs were constructed for each snapshot, and analyzed with our beam-tracer 
algorithm, which is similar to the one described in Kitzbichler & White (2006). We present 
the details of the beam-tracer algorithm in Trenti & Stiavelli (2007). In brief, the pencil 
beam associated with the ACS field of view is traced through the simulation box across 
the different snapshots of increasing redshift. As the uncertainty in the V-dropout redshifts 
is larger than the size of the simulation box, the simulation box needs to be beam-traced 
several times. We carried out 4000 Monte Carlo realizations of the beam-tracing procedure 
by varying the initial position of the beam; the beam was always angled so as to guarantee 
no overlap between different beam-tracing runs. 

Dark matter halos within a beam were flagged, and a simple model •was used to determine 
the number of V-dropouts within a flagged halo. In particular, the intrinsic number of 
galaxies was computed using the halo occupation model description; this gives an average 
occupation number (Wechsler et al. 2001): 

(■ Ngai ) = 9(M - M min ) (1 + (M/Mi) 01 ), (7) 

where 6(x) is the Heaviside step function, M min is a minimum halo mass threshold, M\ 
is the typical scale where multiple galaxies are present within the same halo, and a > 0 is of 
order unity. If M > M min , in addition to the one galaxy at the center of a halo, a suitable 
number of companion galaxies is extracted from a Poisson distribution with mean (M /M{) a . 
A fraction rj of the total galaxies in each halo is finally identified as V-dropouts, so as to 
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take into account that V-dropouts may be detected for only a limited period during their 
evolution. 

The results for different choices of the model parameters are reported in Fig. 8. We 
consider: (i) 77 = 1, Mi — ► +00 (i.e., one galaxy per halo); (ii) 77 = 0.5, Mi — * +00 (i.e., 
one galaxy per halo with 0.5 detection probability); (iii) 77 = 1, M\ = 5 • 10 n M Q /h and 
a = 1 . The value of M rnin is kept as a free parameter and it is adjusted to have the same 
average number of counts in all three cases. Variations of M min are limited in the range 
M min = 2 - 3.5 • 10 w M Q /h. 

The distribution of the number counts probability is very similar in the three cases, 
and the field-to-field standard deviation is a « 25%. This is significantly larger than the 
Poisson uncertainty of sa 7% for a value of ~ 200 counts per field, which is about the number 
of Woe-dropout counts/field that we have detected, after correction for incompleteness and 
selection probability. Depending on the specific model, there are minor differences in the 
computed standard deviation (Act ~ ±2%): for 77 < 1, a smaller M m - m is required; therefore, 
the variance decreases, as lower-mass halos become less clustered. As expected, the largest 
variance results in the case of multiple galaxy occupancy of the halos. 

We also computed the relative variations in expected counts for two fields which are 
separated by about « 10 arcmin, as is the case for the HUDF and the NICP12. field. We 
find that the two fields have essentially independent counts, i.e., the variance of the average 
number of detected V-dropouts decreases by a factor y/2. 

Finally it is interesting to highlight that the narrow-beam geometry of our data, signifi- 
cantly reduces the cosmic scatter with respect to a similarly-sized, spherical volume, which 
has a variance in the number counts of « 50% (see also Somerville et al. 2004). This is due 
to the stretched redshift space of the pencil beam, that probes a variety of environments by 
minimizing the probability of enclosing only a large void or a rich cluster. 


8. Discussion 

We discuss below the implications of our measurement of the faint-end of the z — 5 
LF of star forming galaxies for the evolution of galaxies, in the context of current models of 
galaxy formation. 
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8.1. The z = 3 to 6 Evolution of the Faint-End of the LBGs Luminosity 

Function 

In order to investigate the evolution of the LBG LF over cosmic time, we use the 
published z = 3 and z = 6 LFs by Steidel et al. (1999) and Bouwens et al. (2006), respectively. 
These are plotted in Figure 9. The z ~ 3 LF is K-corrected to 1400 A and adapted to our 
cosmology (M*(z — 3) = —20.9); it is based on identical assumptions as in our fiducial 
LF, with which it can thus be directly compared. The z = 6 LF is based on a ’’cloning” 
algorithm, which is well bracketed by our two different assumptions for the dust extinction 
and the hydrogen intergalactic absorption. 

The possible underdensity of the HUDF, and the associated uncertainty in the faint- 
end slope of the z = 5 LF, is relevant in establishing the evolution of such faint-end through 
the comparison with the z = 3 — 6 estimates. In particular, the upper panel of Figure 9 
shows that the inclusion of the HUDF in the computation of the z — 5 LBGs LF leads to a 
shallower faint-end slope than at z = 3; in contrast, the faint-end slope that we derive using 
the NICP12 field alone for the z = 5 LBGs is very similar to the one measured at that later 
epoch. A similar result is found when the SMC extinction law and the Meiksin intergalactic 
extinction relation are assumed, as shown in the lower panel of the figure. This result is 
thus robust also towards possible variations of the dust extinction properties through cosmic 
time. The possibility remains that the HUDF is, as we discuss above, ~ 30% under-dense 
relative to the NICP12 (and the mean); a renormalization of the HUDF by such an amount 
would actually lead to a slightly steeper faint-end slope at z = 5 than at z = 3. 

Rather similar conclusions are drawn for the redshift evolution of the faint-end of the 
LF in the z = 5 to z — 6 redshift window. Inspecting both panels of Figure 9, our fiducial 
z = 5 faint-end slope is consistently shallower than the z = 6 estimate, for both the HUDF 
and the NICP12 field; the NICP12 field provides however, at z — 5. a similar faint-end 
slope to the one measured at z = 6 when the SMC dust extinction law and the Meiksin 
intergalactic absorption correction are adopted. A similar agreement between the z = 5 and 
z = 6 faint-end slopes would also be. obtained if the HUDF were renormalized to account 
for a ~ 30% underdensity caused by cosmic variance. Note that Bouwens et al. (2006) did 
renormalize the z = 6 HUDF number density, increasing it by about 32%, based on the 
comparison with the GOODS data. 
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8.2. Comparison with Theoretical Predictions 

We compare our observed z — 5 LEG LF with theoretical predictions from the GAL- 
FORM semi-analytical model (Cole et al. 2000), which calculates galaxy formation in the 
framework of the CDM model of structure formation. It includes the assembly of dark matter 
halos by mergers, heating and cooling of gas in halos, star formation from cold gas, feedback 
from supernovae and photoionization, galaxy mergers, and chemical evolution of stars and 
gas. Star formation occurs in two modes: quiescent in disks, and starbursts triggered by 
galaxy mergers. Galaxy luminosities are calculated based on stellar population synthesis, 
and include obscuration by dust, which is calculated self-consistently by radiative transfer 
using the GRASIL model (Silva et al. 1998; Granato et al. 2000). The version of the model 
which we use is that described in Baugh et al. (2005) (see also Lacey et al. 2007a), and 
is based on the A-CDM cosmology, with fi 0 — 0-3, h = Ho/(100k'm,s~ 1 Mpc) = 0.7, and 
a 8 = 0.9. An important difference of this model from the earlier Cole et al. (2000) model 
is the much larger role of star formation bursts at high redshift. In the Baugh et al. (2005) 
model, bursts are responsible for 30% of the total star formation when integrated over all 
redshifts, compared with around 5% in the Cole et al. model, and bursts dominate the total 
star formation rate density at z > 3. This difference in behaviour results from (i) modifying 
the quiescent star formation time-scale in galactic discs, in order to make mergers at high 
redshift more gas-rich; and (ii) allowing triggering of bursts by both minor and major merg- 
ers. The predictions of this model for LBGs at different redshifts will be presented in full in 
Lacey et al. (2007b), here we give only the results for z = 5. The predictions of this same 
model for Lyor-eraitting galaxies at similar redshift, s have already been given by Le Delfion 
et al. (2006). 

The Initial Mass Function (IMF) has a large impact on the predicted z — 5 LF. We 
therefore computed the results for two different assumptions, namely (1) a model where 
all stars are produced with a Kennicutt (1983) IMF, and (2) a model where quiescent 
star formation is described by a Kennicutt IMF, but star formation in the merger-induced 
starbursts occurs with a ’’Top-Heavy* {x = 0) IMF. The Kennicutt IMF has x = 0.4 for 
to < Mq and x — 1.5 above this threshold; it is therefore quite similar to a Salpeter (1955) 
IMF for to > M & . An IMF similar to the Kennicutt IMF has been argued to better fit the 
observational data from our Galaxy and other spiral galaxies than a Salpeter IMF (e.g. Scale 
1998). At constant stellar mass, the far-UV (1000-3000A) luminosity for a Kennicutt IMF 
is within 10-20% of that for a Salpeter IMF. The primary motivation in Baugh et al. (2005) 
for introducing a top-heavy IMF in starbursts was in order to explain the number of sub-min 
galaxies seen at high redshift - with a normal IMF, the predicted number counts of faint 
sub-mm galaxies were too low by more than an order of magnitude. However, a top-heavy 
IMF in starbursts is also motivated by theoretical calculations (e.g. Larson 1998), which 
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suggest that, under the conditions that exist in starbursts, the IMF could be biased towards 
more massive stars relative to quiescent star formation. 

Figure 10 shows the rest- frame 1500A LF at z = 5 predicted by the semi-analytical 
models (Lacey et al. 2007b), compared to our measurements. The error bars on the theoret- 
ical curves show bootstrap estimates of the uncertainty due to the finite number of galaxies 
in the simulated samples; the numerical runs are very CPU intensive, and it is thus currently 
impractical to greatly reduce these error bars by increasing the sample size. The dotted black 
and blue lines in the figure show, for the simple-Kennicutt and Kennicutt+Top-Heavy IMFs, 
respectively, the predicted LFs without dust obscuration. These substantially overpredict 
the number density of z = 5 LBGs over the whole luminosity range sampled by the data, 
but this is not surprising, since the average dust extinctions of star-forming galaxies in the 
far-UV at high redshift are thought to be quite large. 

The solid black and blue lines are the simple-Kennicutt and Kennicutt+Top-Heavy 
IMFs models including the effects of dust extinction. At the bright-end of the z = 5 LF, 
the simple-Kennicutt model with dust extinction substantially underestimate the number 
density of bright (Muv < —19) z = 5 LBGs. The Kennicutt+Top-Heavy IMF model is 
much closer to the observational data. Although this latter model still appears somewhat 
low compared to the observations, this could be explained by relatively small errors in the 
dust extinction, to which the model predictions are quite sensitive. Resolving this issue may 
require better observational estimates of the dust extinctions in LBGs. Another possibility 
would be that either the fraction of mass produced in bursts or the fraction of high mass 
stars in the bursts is even higher than in our Top-Heavy IMF model. 

The Kennicutt+Top-Heavy IMF model also appears to give a somewhat better fit to 
the faint end of the observed z — 5 LF than the simple-Kennicutt model, although the 
difference between the two models at the faint end is only slightly larger than the scatter 
in our estimate due to completeness corrections, input physical parameters, and possibly 
cosmic variance. 


9. Summary and Concluding Remarks 

We have presented the UDF05 project, a follow-up of the HUDF aimed at increasing the 
area of the HUDF and securing two more fields with ultra deep ACS and HI CMOS imaging. 

In this first paper we use the UDF05 to estimate the z ~ 5 Luminosity Function of 
Lyrnan-Break Galaxies down to very faint absolute magnitudes. Specifically, we use the 
NICP12 ACS field of the UDF05, and the ACS HUDF for a comparison, to extend the 
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measurement of the LF of star forming galaxies at z = 5 down to Muoo — —17.1. This 
allows us to constrain the faint-end slope of the z — 5 LF very accurately to a ~ —1.6. 

Attention was paid to correct the raw number counts for selection biases and photometric 
errors, and to estimate the impact on the LF of different assumptions for the properties of 
the z = 5 LEG population. We find that different assumptions about the LBG galaxy 
population and intergalactic absorption result in an uncertainty in the density of LBGs at 
2 = 5 of about 25%. 

Under similar physical assumptions, the HUDF is underdense in z = 5 LBGs by about 
30% with respect to the NICP12 field, a variation which can be accounted for by the expected 
amount of cosmic variance in the pencil-beam volumes that we have probed. 

When added to studies of the faint-end of the LF at lower and higher redshifts, our 
results substantiate a picture in which, once cosmic variance is taken into account, the faint- 
end of the LF of LBGs remains constant through the entire redshift z — 3 to 6 cosmic 
period, during which however the bright Lyman-break galaxy population shows a significant 
evolution. This is in contrast to the results of Iwata et al. (2007), who find only an evolution 
in the faint-end. We note however that, as argued by Ouchi et al. (2004), the selection 
window of Iwata et al. (2007) is rather broad, and might therefore be affected by lower- 
redshift interlopers; this would lead to an overestimation of the bright-end. Also, our results 
are based on data which extends about 3 magnitudes fainter than the Iwata et al. (2007) 
sample. 

The comparison of our data to theoretical predictions suggests that the majority of stars 
of z — 5 LBGs are born with a Top-Heavy IMF in merger-driven starbursts. Still, the bright 
end of the LF is underpredicted by our Top-Heavy IMF model, which could be accounted 
for with a variation in the dust extinction correction. However, at face value, a physical 
interpretation for the discrepancy at the bright end cannot be excluded, i.e. it is possible 
that the fraction of stellar mass produced in starbursts or the fraction of high mass stars in 
the bursts is increasing towards the bright end of the LF. 


Acknowledgements: This work has been partially supported by NASA HST grant 01168. 
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A. Comparison with Published z — 5 LBGs Luminosity Functions 
A.l. The Subaru Deep Field 

Ground-based Subaru observations have been used to derive the z ~ 5 LF of LBGs 
(Iwata et al. 2003; Ouchi et al. 2004; Yoshida et al. 2006). In particular, Yoshida et al. 
(2006) extended the analysis of Ouchi et al. (2004) on the SDF with observations that reached 
0.5 mag fainter and were wider in area by a factor of 1.5. As both authors adopted very 
similar techniques and their number counts of z ~ 5 galaxies are in very good agreement, 
Ave focus on the results by Yoshida et al. (2006) only. Using the combination of Subaru 
Y, R, i\ z' filters, these authors defined two samples of z ~ 5 LBGs, whose reliability was 
tested with spectroscopic data. We can directly compare with their Ri'zf analysis, as this 
filter combination is sensitive to a very similar redshift window {z — 4.5 — 5.4) as our data. 

Following a procedure similar to the one adopted in this paper, Yoshida et al. (2006) 
estimated the selection probability and completeness by adding artificial galaxies to the orig- 
inal images; they did not correct- their samples for photometric errors as we do, but, estimated 
the number of expected low-2 interlopers by simulating galaxies with colors according to a 
photo- 2 catalog from the Hubble Deep Field North (Williams et al. 1996). 

The best-fitting LF derived by Yoshida et al. (2006) is listed in Table 4 (only VHzf 
sample, for the cases of keeping all LF parameters free and of fixing the faint-end slope to 
their z ~ 4 result). In both cases the faint-end slopes -are very steep, i.e., significantly steeper 
than the one that we measure from our 3-magnitudes deeper data. This highlights the need 
for very deep imaging surveys to study the contribution to galaxy assembly at high redshifts 
by galactic systems well below the few bright percent of the star forming population at those 
epochs. 


A.2. A Previous Analysis of the Hubble Ultra Deep Field 

Beckwith et al. (2006) present the first attempt to derive the B435- to Urs-dropout 
LFs in the HUDF. Their approach is very different from ours, as they did not correct for 
observational biases on the selection of dropout sources. Instead, they assumed that these 
biases affect all dropout samples equally, and base on this assumption their study of the 
evolution of the LBGs LF in the z — 4 — 7 redshift window. They extract an '^-selected 
sample from the GOODS vl data to constrain the bright-end of the LF, and derive the 
following best fit parameters: <j>„ = 0.001 Mpc -3 mag -1 , M*(1400 A) = —20.5, a = —1.6 
(fixed). 
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Figure 1 1 shows our fiducial z = 5 LF in comparison with the estimate of Beckwith et 
al. The two LFs are in good agreement, except for one datapoint at M = —18.2. This is a 
very interesting result, as it suggests that the several effects, which we have corrected for in 
our analysis, likely cancel out, thereby validating the Beckwith et al. (2006) assumption, at 
least for the z = 5 LBG population. In a future paper (Oesch et al. 2007, in preparation) 
we discuss whether a similar conclusion can be drawn at z — 6 . 


B. Simulations: Details 
B.l. Idealized-Galaxies Simulations 

Two different types of light profiles were considered for the ideal-galaxies simulations, 
namely an exponential and a deVaucouleur profile. These distributions are expected to more 
or less bracket the whole range of galactic light profiles. Specifically, the surface brightness 
distribution of disk galaxies was approximated by: 


I(x,y) = r— — re 


\/ xr-ja-+y~/b 2 


(Bl) 


2tt a b 

where a and b are the major and minor axes lengths, the inclination angle i is given by 
cos(f) = 7 , and L tot is the extrapolated total galaxy light. The coordinates x and y are 
measured from the centre of the galaxy. The disk was assumed to be infinitely thin, which 
would lead to unrealistic shapes for edge-on systems. To avoid this unphysical situation, 
galaxy shapes were simulated with inclinations up to a maximum of 80°. To retain a flat 
distribution' of inclinations in the entire 0° to 90° range, inclinations > 80°, drawn from an 
assumed fiat distribution, were set equal to 80°. 


The deVaucouleur profile adopted to simulate spheroidal- type galaxies reads: 


m = L ^ i ^? e ' 7 ' 6692(r/re//)I/4 ( B2 ) 

with r the radial distance as measured from the centre of the galaxy, and L tot the 
extrapolated total galaxy light. Spherical symmetry was assumed for the deVaucoleurs 
galaxies. The r 1//4 -profile is very steep at small radii, and thus these types of galaxies are 
detected down to a fainter limiting magnitude with respect to disk galaxies. The total light 
is however distributed over a more extended region- and therefore, the amount of flux lost 
in the low-surface brightness wings is larger - than for an exponential profile. 
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The idealized galaxies were convolved with the ACS PSF before being added to the 
NJCP12 and HUDF images. The PSF was derived from a relatively bright but unsaturated 
star in the field. Tests with a PSF produced with TinyTim, downgraded to the measured 
FWHM of the observed PSF, were also performed, and produced similar results. Poisson 
noise from the galaxy counts was not added to the simulated galaxies, since at the faint 
magnitudes of this study the sky background is the dominant source of Poisson noise. Tests 
were nonetheless performed, which confirmed that the galaxy shot noise had no impact on 
the results. 


B.2. Dimmed-Galaxies Simulations 

The idealized-galaxies simulations allow us to relate the intrinsic properties of galaxies 
to how these would be observed with our instrumental setup. They suffer however from 
the limitation of assuming perfectly smooth light distributions, as described by the adopted 
analytical descriptions. This is clearly not the case for real galaxies at any epoch, and 
certainly not for the UV light emerging from the star forming z ~ 5 galaxy population. 

In order to have a benchmark to gauge the incompleteness in our samples on the basis 
of more realistic assumptions concerning the galaxy light distributions, we performed an 
additional set of simulations, entirely based on the galaxies that were detected in the NICP12 
and HUDF fields. Specifically, we dimmed the ACS images by specific amounts, and added 
these dimmed images to the corresponding original frames. Offsets of a few pixels were 
introduced before co- adding the dimmed and original images. In order to maintain the S/N 
ratio properties of the original data, the co-added images were appropriately rescaled. More 
specifically: 

1. A copy of the original image was divided by a factor /; 

2. This /-scaled version of the data was shifted by a few pixels relative to the original 
parent frame and added to the latter; 

3. The resulting composite image was di Added by + ljff, to retain the noise prop- 
erties of the original image; 

4. The rms map of this final ” dimmed image” was obtained from the rms map of the 
original image by computing: 

\] v'thifj 'f + tU 
x/1 TTJp 


rms 
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The procedure leads to an effective dimming factor for the galaxies’ fluxes of n^i m = 
y/\ -f ,/’ 2 . We performed four separate experiments with four different shifts and dimming 
factors ridim — \/2, 2,4, 6. We used these ’’dimmed-galaxies” simulations to quantify, in an in- 
dependent and complementary way relative to the ” idealized- galaxies” simulations discussed 
above, the completeness of our catalogs, as well as the scatter in measured magnitudes and 
colors, as a function of galaxy luminosity. 
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Fig. 1. — The planned UDF05 fields relative to the original HUDF ACS and NICM OS- 
parallel pointings. The larger and smaller squares correspond to the ACS and NICMOS 
pointings, respectively. The V-dropout study that we present in this paper is based on the 
UDF05 ACS observations covering the NICP12 field, and the HUDF. Due to the two-gyro 
operation of the HST, the UDF05 adds a new deep-imaging NICMOS field to the original 
HUDF pointings (NICP5) . 
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Fig. 2. — Sjmthetic SED of a starbursting galaxy at z — 5 with z$ 50 = 25.5 mag. The optical 
passbands used for the UDF05 program are overplotted. The vertical dashed lines indicate 
the Ly a edge at redshifts 4.5, 5 and 5.5, respectively. 
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Fig. 3. — Evolutionary tracks of different types of galaxies, and location of the stellar 
sequence, in the Vqqq — ^750 versus 1775 — zsso color-color diagram. The tracks are con- 
structed using the Bruzual & Chariot (2003) population synthesis models. The thick black 
line corresponds to the adopted selection (Eq. 1 - 3). The thin solid lines show the 
tracks of 200 Myr old, continuously star-forming galaxies at z = 4 with E(B — V) — 
0 (blue), 0.15 (green), 0.3 (red) for a Calzetti et al. (2000) dust extinction relation. Lower 
reclshift galaxy types are plotted with dotted, dash-dotted and dashed lines (z — 0 — 4), cor- 
responding to local irregulars, Sbc spirals, and a single stellar population with an age of 3.5 
Gyr, respectively. The latter was chosen because it corresponds to a maximally old stellar 
population at 2 ~ 1.5, the redshift at which single-burst early-type galaxies contaminate the 
selection. The small crosses correspond to the 131 Galactic stars from the Pickles (1998) 
library. 
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Fig. 4. — The predicted redshift distribution for V^g-dropout -dropouts as derived assuming 
a non-evolving LF and synthetic SEDs (see text). The total interloper fraction is estimated to 
be 23% and is primarily contributed by very reddened galaxies with either strong emission 
lines or an old stellar population component. The inset panel zooms over the primary 
redshift range of interest. The non-uniform distribution of the histogram highlights the fact 
that different SEDs are responsible for the signal in different redshift intervals. 
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Fig. 5. — Estimated detection completeness over the range of interest in the two fields as a 
function of observed magnitude from the idealized-galaxies (solid line) and dimmed-galaxies 
(dashed line) simulations. The completeness does not reach 100% even at the brightest 
magnitudes; this is due to blending of sources bj r the SExtractor algorithm. 
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Fig. 6. — Top — The NICP12 selection probabilities for the different assumptions on the 
SED distributions and intergalactic hydrogen absorption discussed in the text, for the case 
of a bright (25 — 25.5, blue lines) and faint (28 — 28.5, red lines) magnitude bin. The solid 
(dot-dashed) and dashed (dotted) lines refer to the starburst (SMC) dust extinction relation, 
assuming the intergalactic hydrogen absorption prescription of Madau (1995) and Meiksin 
(2006), respectively. Bottom — The same for the HUDF. 




Fig. 7. — Total V e g LF obtained combining the NICP12 and HUDF data (this work) with 
the SDF LF (Yoshida et al. 2006, i&V-sample). Red squares and blue circles correspond to 
the NICP12 and HUDF V e g LFs respectively. Green crosses are the V^-LF by Yoshida et 
al. (2006). The errorbars indicate Poissonian errors only. The solid line shows the best- fit 
Schechter function to the total LF (see text and Table 4). 
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Fig. 8. — Probability distribution of the V-dropouts number counts in one ACS field (top) 
and in the two HUDF + NICPIC12 fields (bottom), accordingly to different halo occupation 
distribution models. Black line: one galaxy per halo, q = 0.5; red line: one galaxy per halo 
q — 1; blue line: HOD model with Mi = 5 • 10 n M®/h and q = 1. 
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Fig. 9. — Evolution of the rest-frame 1400 A LEG LF in the redshift range z ~ 3 — 6. 
Top panel: The LFs assume the starburst dust extinction relation and the intergalactic 
absorption prescription of Madau (.1995). Our fiducial .z = 5 LF is represented by the black 
circles (NICP12 data) and the black squares (HUDF data), plus the black crosses, which 
represent the z = 5 SDF LF of Yoshida et al. (2006). In blue we show the z ~ 6 LF of 
(Bouwens et al. 2006), and in red the z ~ 3 LF of (Steidel et al. 1999). Solid and dotted 
lines are used to indicate the magnitude baseline over which the LF were measured and 
extrapolated, respectively. Bottom panel: A similar plot as above, in the assumption of 
a SMC dust extinction la,w and a Meiksin (2006) intergalactic absorption correction. The 
parameters for the LFs are listed in Table 4. 
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Fig. 10. — Comparison of our fiducial 2 = 5 LF (red solid line: fit to the NICP12. HUDF 
and SDF data; red circles: NICP12 data; red diamonds: HUDF data; black crosses: SDF 
data) with the prediction of Lacey et al. (2007). The simple- Kennicutt and Kennicutt+ : Top- 
Heavy IMF models are shown in black and blue, respectively; dotted and solid lines indicate 
unobscured and dust-extincted models. 
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Fig. 11. — Comparison of our LF to the published HUDF LF by Beckwith et al. (2006). Our 
best fit to the combination of the HUDF data (red diamonds), NICP12 data (black squares), 
and the SDF data (blue crosses) is shown as a black line, and the estimates of Beckwith et al. 
are shown in green (symbols and line for the measurements and best fit Schechter function, 

respectively). 




41 


Table 1. Properties of the UDF05 and HUDF Observations 


Field 

Filter 

HST Orbits 

Exposure Time [s] 

10cr-magnitude a [AB] 

NICP12 

Deo6 

9 

2U600 

28.49 

NICP12 

'irrs 

23 

54’000 

28.44 

NICP12 

-*850 

69 b 

168’000 

28.47 

HUDF 


56 

134’900 

29.38 

HUDF 

Deoe 

56 

135’'300 

29.81 

HUDF 

^775 

144 

347T00 

'29.43 

HUDF 

-2850 

144 

346 ’600 

28.73 


a within an aperture of O'/I 5 radius 
b 70 planned; one lost due to loss of lock 
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Table 2. Number of detected LBGs and corrections in the two fields 


'^850 t 

N catalog 

bootstrap 

fint N corrected l 

^-'Idealized 

(-'Dimmed 


NICP12 

25.25 

6 

4 

0 

4.0 

0.95 

0.95 

25.75 

4 

6 

0.12 

5.3 

0.93 

0.93 

26.25 

12 

12 

0.1 

10.8 

0.92 

0.92 

26.75 

14 

15 

0.14 

12.9 

0.91 

0.90 

27.25 

12 

22 

0.19 

17.8 

0.88 

0.85 

27.75 

29 

28 

0.24 

21.3 

0.76 

0.79 

28.25 

24 

32 

0.28 

23.0 

0.64 

0.68 

HUDF 

25.5 

3 

3 

0 

O 

O 

0.92 

0.92 

26 

2 

4 

0 

4 

0.91 

0.91 

26.5 

6 

9 

0.04 

8.6 

0.90 

0.89 

27 

23 

19 

0.09 

17.3 

0.87 

0.88 

27.5 

18 

21 

0.13 

18.3 

0.84 • 

0.86 

28 

24 

28 

0.18 

23.0 

0.78 

0.81 

28.5 

34 

35 

0.10 

31.5 

0.68 

0.70 

29 

23 

27 

0.08 

24.8 

0.47 

0.53 


(The columns list, respectively, msso- central bin magni- 
tude; N catalog- uncorrected number of color-selected LBG sources; 
N, bootstrap : sources after bootstrapping; fi nt : expected interloper 
fraction due to photometric errors; N correcte d' final corrected Veo6- 
dropout sources; Oid ea x Vie d- completeness derived from the idealized 
galaxies simulations; Coimmecb completeness derived from the dim- 
ming simulations 





Table 3. Impact of Physical Parameters on the LF (Eq. 4; Constant M* = —20.7)1 
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Table 4. z ~ 5 Luminosity Functions 
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